Self-Starting Micromotors in a Bacterial Bath 
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Micromotors pushed by biological entities, like motile bacteria, constitute a fascinating way to 
convert chemical energy into mechanical work at the micrometer scale. Here we show, by using 
numerical simulations, that a properly designed asymmetric object can be spontaneously set into 
the desired motion when immersed in a chaotic bacterial bath. Our findings open the way to 
conceive new hybrid microdevices exploiting the mechanical power production of bacterial organisms. 
Moreover, the system provides an example of how, in contrast with equilibrium thermal baths, the 
irreversible chaotic motion of active particles can be rectified by asymmetric environments. 



Ensembles of animate organisms behave in a very rich 
and surprising way if compared to inanimate objects, as 
atoms or molecules in a gas or a liquid. Everyone has 
been amazed by the cooperative motion of birds in a 
flock, fishes in a school or wildebeest in a herd p], Q. 
Also at the micrometer scale elementary living organ- 
isms, like bacterial cells, show an extraordinary variety of 
behaviors, such as collective motions 0,0,0,0], complex 
chemical- mediated motility or chemotaxis [7[ , spatiotem- 
poral patterns self-organized structures [§], biofilms 
formation [loj]. An important peculiarity of animate or- 
ganisms is the fact that they can be self-propelled, using 
a variety of different mechanisms for this purpose [11]. 



Motile cilia and turned flagella are two example of evo- 
lutionary tricks adopted by living organisms to accom- 
plish the hard task of swimming at low Reynolds number 
UJ. One can think about such ensembles of organisms 
as open systems, with a net incoming flux of energy (pro- 
vided by nutrients) stored and converted into mechani- 
cal motion by irreversible processes happening inside the 
cell body. The resulting dynamics breaks time inversion 
symmetry so that asymmetric environments can result in 
directed motions which, in equilibrated Hamiltonian sys- 
tems, would be forbidden by detailed balance [H, EH- A 
natural question then arises: is it possible to rectify such 
a non equilibrium dynamics to propel microdevices? 

Biological molecular motors constitute a fascinating 
mechanism to generate motion at the nanoscale [HI, [l6| . 
When larger, micron sized, structures need propulsion 
the preassembled motor units found in unicellular motile 
organism may offer several advanta ges over isolated pro- 
teins. In a recent experiment [~L7l Il8| bacterial driven 
micromotors have been assembled by biochemically at- 
taching motile bacteria to a microrotary motor. Such 
procedures require the construction of narrow tracks to 
induce a unidirectional binding of bacterial cells on to the 
moving rotor with a consequent increased complexity in 
designs and limited number of working bacteria. 

Here we numerically show that a properly designed 
asymmetric motor immersed in a chaotic bacterial bath 
can be spontaneously set into the desired motion. Our 
numerical findings suggest the possibility to construct 



new opportunely shaped microdevices able to exploit the 
propelling power of motile bacteria. 

Spinning a bundle of helical flagella, bacteria like E. 
coli, may swim along their body axis with speeds of or- 
der 10 body lengths per second jTj|. Decorrelation of 
velocity may occur via four different mechanisms: tum- 
bling, Brownian motion, mechanical interactions and hy- 
drodynamic interactions. The first mechanism is a spon- 
taneous tumble produced by a temporary reversal in the 
spinning direction of the flagellar motor [20]. Brown- 
ian motion can also be effective in producing diffusion 
of orientation and hence of propelling direction. Inter- 
actions with other bacteria can be mechanical, by direct 
contact, or hydrodynamic, via flow currents produced by 
the swimming motions. Trying to mimic the behavior of 
an elongated E. coli cell with a minimal model we only 
retain the two most effective mechanisms that are tum- 
bling and mechanical interactions. Hydrodynamic inter- 
actions, occurring only through dipole or higher order 
multipoles, turn out to be effective only over short dis- 
tances where mechanical interactions between elongated 
bodies are much more effective in reorienting the bacte- 
ria. We directly checked that including hydrodynamic 
interactions has a negligible effect on the mean squared 
displacements and on its crossover from ballistic to dif- 
fusive regimes. 

Each cell is represented by an instantaneous position 
and an orientation e$ pointing in the free swimming 
direction. The elongated hard body of the cells (length I 
and thickness a) is modeled by the sum of p short range 
repulsive potentials centered at equally spaced locations 
along the cell axis rf = + d^ei with /3 = l,p and 
dP = (Z - a) (2/3 -p - l)/(2p - 2). The neighboring cells 
will then act on the ith cell with a system of forces 
applied at rf : 
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FIG. 1: (Color online) Rotary micromotor in a bacterial bath. Snapshots at three different simulation-times, t = 10, 12 and 
14 s. Each bacterium is represented by a spherocylinder (with aspect ratio 1/2) whit a white head pointing in the direction 
of the self-propelling force. The arrow at the center of the gear evidences the counterclockwise rotation at an average angular 
velocities coo ~ 0.21 rad/s. 



To such intercellular forces we added intracellular forces 
consisting of a constant linear propelling force fo (direct- 
ing along &i) which is only active in the running state 
and a random torque T r which switches on during the 
tumbling state. The probability per unit time to switch 
in a tumbling state is constant and such as to give an 
average free run length of 10 cell lengths [l9[. Introduc- 
ing the state variable 8i which is in the running state 
and 1 during a tumbling event, the net forces and torques 
acting on the ith cell read: 



/ e,(l-^) + E F; 
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For the subsequent motion the rigid cell body is modeled 
as a prolate spheroid of aspect ratio a = a/ 1. Therefore 
the center of mass and the angular velocity are [2l|: 



Vi = Mi • Fj 
= K r T, : 



(5) 
(6) 



where 



We choose I as the unit length, r = l/vo as the unit of 
time (where = m \\fo is the free swimming velocity) 
and m\\ as unit of mobility. When not specified physical 
quantities will be expressed in reduced units. A planar 
geometry will be investigated in a box LxL with periodic 
boundary conditions. We will specialize to the case of 
TV = 1092 bacteria with number density p = N/L 2 = 
0.945, aspect ratio a = 1/2 and potential parameters 
p = 2, n = 12. Mobility values are m\\ = 1, m± = 0.87, 
k± = 4.8 (fc|| does not enter in the equation of motion, 
as Ti is perpendicular to in the planar geometry). We 
consider a micromotor immersed in the bacterial bath. 
The asymmetric micromotor is a gear with a sawtooth 
profile whose center of mass is kept fixed at the center 
of the box. The motor is free to rotate around its axis. 
Each of the p force centers, describing a single bacterium 
body, interacts with boundary walls through a force of 
the form in Eq. [2j where r is a vector perpendicular to 
the wall connecting the p-centers to a point located at a 
distance a/2 behind the wall. 

The resulting cell-boundary forces produce further 
force and torque terms in Eq.s[5]and[6j and a net fluctu- 
ating torque on the gear motor, whose angular velocity 
is then 



= raii^e* + ra_L (1 - (7) 
K* = feline* + k± (1 - (8) 

We choose the force coefficient A in such a way that two 
bacteria facing head to head on the same line would be 
in equilibrium at a distance a = al: 

A/a n+1 =f ^A~ / a n+1 (9) 



Slg = Kg Tg (10) 

where T g is the torque exerted by bacteria on the gear 
whose rotational mobility is K g . We consider a gear with 
8 teeth and internal (external) radius Ri nt = 5 ( R ex t = 
8). The gear mobility is estimated as that of a disk [2l( of 
radius 6.5: K g = 1.9 • 10 -3 . Equations of motion ([5]), (|6|) 
and ([TP]) where numerically integrated by Runge-Kutta 
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FIG. 2: (Color online) Angular velocity uj (in rad/s) of the mi- 
cromotor as a function of time (in seconds). Black- line refers 
to a single run, red-line is the average over 100 independent 
runs. After a short transient regime (due to initial configura- 
tion of bacteria) , a fluctuating velocity around a mean value 
coo ~ 0.21 rad/s is observed. Inset: same as main plot for a 
symmetrically shaped micromotor. 



method [22[ for 2 • 10 5 steps (with time step St = 10" 3 ). 
At t = the bacteria are uniformly distributed in the 
space outside the external disc or radius r ext . 

We found that the micromotor starts to move spon- 
taneously under the effects of pushing bacteria. A net 
unidirectional motion is observed, with a fluctuating an- 
gular velocity around a non-zero mean value. In Fig. 1 
we show snapshots of the bacterial bath with a rotary 
micromotor at three different times, t = 10, 12 and 14 s 
(physical units are obtained considering realistic values 
I = 3 /im, = 30 jitm/s). It is evident a densification 
process close to the device's boundary, in agreement with 
recent studies on self-propelled cells in confined environ- 
As a result a net rotary counterclockwise 
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ments 

motion of the gear during time takes place. The instan- 
taneous angular velocity uo of the motor as a function of 
time is shown in Fig. 2, where black-line refers to single 
run, while red-line is the average over 100 independent 
runs. After a short transient clockwise rotation the sys- 
tem reaches a stationary regime with a fluctuating pos- 
itive (counterclockwise) angular velocity around a non- 
zero average value c^o — 0.21 rad/s, corresponding to 2.0 
rpm . The total torque on the device can be estimated 
around 17 pN /im (assuming m\\ = 59 /im/pN s). The on- 
set of a directed rotation can be understood by analyzing 
the collision of a single bacterium with the rotor bound- 
ary. When a bacterium touches a rotor edge it will exert 
a force given by the projection of the propelling force onto 
the surface normal (black arrows in Fig. 3). The same 
force will act on the bacterium producing a net torque 
that will align the cell body along the edge. Depending 
on the sign of the incident angle measured from the wall 
normal, the bacterium will then quickly leave the gear 




FIG. 3: (Color online) Sketch of the collision of a single bac- 
terium with the rotor boundary. Arrows are the forces ex- 
erting by the bacterium onto the rotor, a) Bacteria coming 
from the left area with respect to the normal leave the gear, 
b) Bacteria from the right get stuck at the corner exerting a 
torque on the rotor. 



back into the solution (Fig. 3a) or get stuck at the corner 
exerting a torque on the rotor (Fig. 3b). The same rea- 
soning applies for both the long and short edges. Most of 
the collisions, however, will occur on the long edge con- 
tributing a transient negative torque which explains the 
short times negative dip in rotor angular velocities. It is 
worth noting that the elongated form of the bacteria is 
not essential for the observed effect, as the same directed 
motion of the micromotor also occurs in the presence of 
"spherical" bacteria, i.e. with aspect ratio a = 1. The 
shape of the motor, instead, plays a crucial role. Indeed, 
simulations performed with a symmetric gear (with sym- 
metrically shaped teeth) produce on average an immobile 
motor, whose angular velocity fluctuates around zero (in- 
set in Fig. 2). The asymmetry is then a basic ingredient, 
as observed in many other thermal ratchet mechanisms 
discussed so far in the literature 
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Given the above mechanism, one expects that the 
torque T g exerted by bacteria would increase as the 
square of the size R of the rotor as both perimeter (and 
hence applied forces) and moment arm increase linearly. 
On the other hand, the rotational mobility of the gear 
Kg decreases as 1 / R 3 resulting in an average angular ve- 
locity decreasing as 1/R. The maximal work that can 
be extracted from the bath is obtained when an exter- 
nal reversible system applies an opposing torque equal 
to Tg/2. The extracted mechanical power is then given 
by TgKg/A and increases with R. Therefore in a planar 
geometry, a 2D array of small gears would perform bet- 
ter, in terms of usable power, than a single big one. The 
dependence on bacterial concentration is also non trivial 
due to the interbacterial interactions that could result 
in reduced motilities at high packing fractions. We note 
that, however, the observed directed motion of the rotor 
is a quite robust effect with respect to the variation of 
different physical parameters, as the density of bacteria, 
their aspect ratio, the shape of the asymmetric rotor, its 



size, the boundary conditions: a quantitative discussion 
on the role of different parameters will appear in a forth- 
coming paper. 

Our main point here is to demonstrate that, in contrast 
to thermal baths of passive particles, useful work can be 
extracted from the chaotic motion of a non-equilibrium 
suspension of active objects. This behavior reminds of 
the ratchet effect or Brownian motors [26] , in which out- 
of-equilibrium systems undergo a rectification process in 
the presence of some asymmetric potential or device. 
More specifically, in an equilibrated Hamiltonian system, 
there's no entropy production and time reversal symme- 
try guarantees that any trajectory has the same probabil- 
ity than its time reversed, so that no systematic directed 
motion can be observed on average. On the other hand 
when a self propelled particle collides to another (or to a 
boundary), the forces they exchange is not just the repul- 
sion of their rigid bodies, but there are also the forces gen- 
erated by the propelling units. Such forces are directed 
along the incoming directions of the two particles and 
therefore would change sign upon time reversal, while 
particles repulsion wouldn't. Time reversed trajectories 
are then incompatible with the assumed dynamical laws. 
From a thermodynamic viewpoint such irreversible dy- 
namics reflects the constant entropy production involved 
in the chemico-physical processes driving the propelling 
unit, such as the flagellar rotary motor of E-coli. Once 
time inversion symmetry does not hold a spontaneous 
directed motion is allowed whenever a spatial inversion 
symmetry is broken [2^ ]. 

In conclusion, we have shown that it is possible to con- 
ceive opportunely shaped microdevices that can move in 
a directional way when immersed in a bath of motile mi- 
croorganisms. In particular, we numerically show that 
a rotary micromotor, consisting of an asymmetric gear 
in a bath of E. Coli bacteria, spontaneously sets into 
a unidirectional rotational motion at an average speed 
of a few rpm. Using asymmetrically shaped boundaries 
also linear translatory motions could be obtained and 
bacterial driven transport could be achieved by self as- 
sembly of bacteria along the particle's boundary. Re- 
markably, when coupled to an external reversible device, 
a net amount of useful energy could be extracted from 
the chaotic motion of a bacterial bath. Our findings can 
open the way to new and fascinating applications in the 
field of hybrid bio- microdevices engineering, and also pro- 
vide new insight in the more fundamental aspects of non- 
equilibrium dynamics of active matter. 
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